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Abstract 

We develop a formalism for the calculation of the ground state energy 
of a spinor field in the background of a cylindrically symmetric magnetic 
field. The energy is expressed in terms of the Jost function of the associated 
scattering problem. Uniform asymptotic expansions needed are obtained 
from the Lippmann-Schwinger equation. The general results derived are 
applied to the background of a finite radius flux tube with a homogeneous 
magnetic field inside and the ground state energy is calculated numerically 
as a function of the radius and the flux. It turns out to be negative, remain- 
ing smaller by a factor of a than the classical energy of the background 
except for very small values of the radius which are outside the range of 
applicability of QED. 
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1 Introduction 



The ground state energy of the spinor field in the background of a magnetic field 
had been investigated since the early days of Quantum Electrodynamics (QED). 
So, for instance, the effective potential in a strong magnetic field is well known. 
For a weak field it takes increasing positive values, for stronger fields it turns 
down and for B — > oo it becomes 



12tt 2 m e K J 

where m e is the electron mass. For details see e.g. HJ. 

Actual interest in this topic results, for example, from the symmetry restora- 
tion due to a magnetic field in the electroweak theory or from the influence of 
that field on the character of the phase transition 0. While most work had 
been done in homogeneous magnetic fields, mainly because in that case explicit 
formulas exist, the extension to inhomogeneous fields is of interest. 

Some work is this direction has already been done. For example for a flux 
where the magnetic field is concentrated on the surface of the tube, the fermion 
determinant was calculated in ||. However, in that case the classical energy 
is infinite. In |4[], in (2+1) dimensions a magnetic field homogeneous in the y- 
direction and with a special shape in the x-direction allowing for explicit formulas 
had been considered. The result for the ground state energy per unit area of the 
planes was expressed in quite elementary functions, which allowed the discussion 
of the total energy within the family of fields considered. It was shown, that 
the system is driven towards a uniform magnetic field. An extension of these 
formulas to the (3+1) dimensional case was given in |JJ. 

Furthermore, there exists a number of investigations of the density per unit 
volume e(r) of the ground state energy, where r is the radial coordinate in the 
perpendicular plane, in the background of an infinitely thin magnetic flux tube. 
Partly, this is motivated by the close relation to the Aharonov-Bohm effect. The 
first investigation of this kind was done in || for the ground state energy density 
of a scalar field, later reconsidered and generalised to the spinor case in |7|], 
see also ||. Extensively investigated is the (2+l)-dimensional case, see || and 
papers cited therein. Also, there are similar investigations in the background of 



an infinitely thin cosmic string [10|. The calculations in the background of an 
infinitely thin magnetic flux tube have the drawback that the energy density per 
unit volume cannot be integrated to get the energy density £ per unit length due 
to the singular behaviour near the string, e(r) ~ r~ 4 , which follows already from 
dimensional reasons. 

In addition one might consider the combined effect of boundaries and back- 
ground fields. This has been started in where imposing spectral boundary 



conditions the spinor field was considered in a finite region of space in the back- 
ground of an Aharonov-Bohm flux string. 



1 



In the present paper we continue the consideration of inhomogeneous mag- 
netic background fields and calculate the ground state energy of the spinor field 
in the background of a straight magnetic flux tube of finite radius R, more ex- 
actly, the energy density per unit length. The reason to consider a flux tube of 
finite radius R is that the associated classical energy is finite and the dependence 
of the total energy when R varies while the flux is fixed can be analysed. The 
interesting question in this context is if some radius R m exists where the complete 
energy, i.e., the sum of classical enegy of the magnetic field and the ground state 
energy of the spinor field, is minimized and the magnetic string gets stable. In 



doing this analysis we use and generalise the formalism developed in [12] and [13 



for a smooth scalar background field. In section 2 we start discussing in detail 
the renormalization of the ground state energy. We will normalize the energy in 
such a way that it vanishes for the electron mass m e — > oo, the resulting massless 
limit will also be considered. The needed counterterms and the subtraction em- 
ployed is then elegantly described using the heat-kernel language. After having 
explained in detail the renormalization procedure we express the ground state 
energy in terms of the Jost function of the associated scattering problem. The 
procedure developed in [12, 13] consists of adding and subtracting the uniform 
asymptotic expansion of the Jost function. This is done in section 4 using a 
perturbative expansion of the Lippmann-Schwinger equation. Various details of 
this calculation are relegated to the Appendices A and B. In the remaining part 
of the paper, section 5 and 6, an analytical as well as numerical description of 
the ground state energy for spinors in the presence of a finite radius flux tube 
is performed. Some details can be found in Appendix C. In the conclusion we 
summarize the main results of the paper. 



2 Basic formulas and the renormalization 

The considered background is a straight magnetic flux tube of finite radius R, 
i.e., the magnetic field 

B(x) = ±- h(r) e z , (2) 
Zn 

where h(r) is a profile function with compact support in the radial variable 
r = y/x' 2 + y' 2 in the plane perpendicular to the tube. By the normalization 
Jq 00 drrh(r) = 1,0 has the meaning of the flux inside the tube. The correspond- 
ing vector potential can be chosen to be 

m - I 2fi2 e - . (3) 

The profile functions h(r) and a(r) are connected by the relation h(r) = a'(r)/r. 
Below, we will use the special case of a homogenoeus magnetic field inside the 
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tube, where these functions read 

h(r) = ^-Q(R-r), a(r) = ^Q(R - r) + 6(r - R). (4) 
R z R z 

In this case the solutions of the field equations can be expressed in terms of Bessel 
and hypergeometric functions. 

Another choice could be a magnetic field concentrated on the surface of the 
cylinder, h(r) = 8(r — R)/R, where the solutions can solely be expressed in terms 
of Bessel functions. But in that case the classical energy of the background is 
infinite. 

The classical energy of the background (per unit length of the string) is 



C class __ 

t "2 



±Jdx B 2 = ^-Jdrr h(rf. (5) 



o 

The ground state energy of the spinor field in that background is given by 

£ = ~ E W" 2s - (6) 

(n,e) 

Here, the minus sign in front of the rhs. accounts for the spinor obeying an- 
ticommutation relations, s (s > 2 and s — > in the end) is the regularisation 
parameter in the zetafunctional regularisation which we use here and ji is the 
arbitrary dimensional parameter entering this regularisation. In fact, due to the 
translational invariance along the axis of the flux tube, we have to consider the 
energy density per unit length. This will be taken into account below. 
The e( n ,e) are the eigenvalues of the Hamiltonian 

U = -i 7 V (d/dx j - ieAj(x)) + 7°m e (7) 

which follows from the Dirac equation. Here, e = ±1 is the sign of the one par- 
ticle energies e( n ,e) f° r the particle respectively the antiparticle which themselves 
are chosen to be positive. All other quantum numbers are included into (n). 
Furthermore, the index j denotes spatial indices only and summation over it is 
included. 

For the renormalisation we follow the standard procedure using the heat kernel 
expansion. The ground state energy can be expressed by the keat kernel K(t) of 

n\ 

L 2s J dt t s-3/2 



£= ~J T(rw) m 

with the asymptotic expansion for t — > 



„— tm 
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where d is the dimension of the manifold, d = 3 in our case. 

The heat-kernel coefficients a n for the operator under consideration are well 
known. The relevant operator 7i 2 reads explicitly, 

H 2 = + -a ij Fij + m 2 , 



with = | [7^7^] and the leading coefficients can be found for example in [14 

For n = we note that the coefficient is independent on the background field 
and corrresponds to the contribution of the empty Minkowski space. We drop 
this contribution without further comment. The coefficient a\ is zero and for a 2 
the general formula reads 



12 J 



a 2 = Tr I dx(-^F 2 + -(a^F tj ) 2 



Here, the trace is over the spinors. The integration along the axis of the flux tube 
gives the corresponding volume by which we have to divide. So the following 
formulas have to be understood always as densities with respect to this axis. 
The trace in ([8]) can be carried out and by means of (|3]) we arrive at 

oo 

a 2 = ^5 2 [drr h{r) 2 , (9) 




MS (io) 



where the notation 



is introduced. Here, e is the electron charge, and we can rewrite this relation as 
5 = y/acj)/ (2tt) , where a is the fine structure constant. 

Using the heat kernel expansion it can be shown that the divergent part of the 
ground state energy results from the contribution of the heat kernel coefficients 
a n with n < 2. We define 

32n 2 \s ml) 
Then the renormalised ground state energy is given by 

Here the limit s — > can be performed because the pole part is subtracted. In 
general, the definition of £ dlv is not unique. By the definition (|11|), the normali- 
sation condition 

£ rcn ->0 for m e ^oo (13) 
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is assured. This normalisation condition is natural as it implies that a very 
massive field should not show quantum fluctuations. On the other hand, it fixes 
the arbitrariness which came in with the parameter /i in (|6|). 

Along with the subtraction of £ dlv from £ it must be added to £ class . This is 
equivalent to a renormalisation of the flux according to 

(14) 

From this renormalisation procedure it is possible to determine the leading 
asymptotic behaviour of the renormalised ground state energy when the radius 
R of the flux tube tends to zero. In fact, this could have been done using the 
arguments given in |T5[ where the general scaling behaviour of the Casimir energy 




was investigated. The point is simply that the regularized ground state energy (y) 
has a series expansion with respect to powers of the mass m e . Note that this is not 
affected by the zero modes, which are present here (see |1|), because we consider 
the ground state energy and not the determinant of the operator V. Then, by 
means of (0) and ( |TT]) we subtract a contribution containing lnm e . Therefor, 
the renormalized ground state energy £ ien becomes for m e — > proportional to 
lnm fi . Now, for dimensional reasons it can be written as 



gren ^ f ' ( Rm c 



R 2 

where / is some function of the dimensionless combination Rm e . Consequently, 
for R — > the behavior must be 



lQir 2 Rm e 
and, for instance, with the background ([|) 

This behavior will be confirmed below in the course of the explicit calculations. 
In fact, this behavior follows already from the heat kernel expansion and has 
the same origin as the asymptotics (|l|). For the same reason it can obviously be 
improved using the renormalisation group. Because the theory is infrared free 
and the limit R — > corresponds to high momenta it is impossible to make R 
too small. 

When comparing ([15]) with the classical energy of this background, 

e*- = ^. us) 

one could think that the complete energy can be made negative for sufficiently 
small R. However, this would require Rm e < exp(— 6n 2 /a), which is far outside 
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the range of applicability of QED and also ruled out by the renormalization group 
argument. 

Note that in the massless case the renormalization scheme must be modified. 
As discussed in [19|, for a 2 7^ which is the case here generically (a 2 = means a 
vanishing background, cf. @) and which represents the conformal anomaly, there 
is no transition for the renormalised ground state energy £ ren from m c 7^ to 
m e = 0. Usually, in the massless case in the definition of S dlv like ([0]) instead of 
the mass, another parameter, the renormalization scale A, is used. Therefor the 
renormalised ground state energy contains a nonuniqueness proportional to the 
heat kernel coefficient a 2 , i.e., proportional to the classical energy (see formulas 
d) and ©). 



3 Ground state energy expressed in terms of the 
Jost function 

We express the regularised ground state energy (|j) in terms of the Jost function 
of the scattering problem associated with the operator H, Eq. (0). Because the 
background is translationally invariant along the third axis, by means of 



we rewrite the corresponding Dirac equation using the chiral representation of 
the gamma matrices in the form 



p + L - m c cr 3 p 3 a 3 \ /$ 

P3°3 p + L + m c a 3 J \ip 



0, (17) 



with L = iY^^i^iid/dx 1 + ieAi). It is sufficient to consider p 3 = and one of 
the two decoupled equations. Then by means of 



igi(r) e - i(m+1) ^ 



g 2 (r) e~ im * 
(m = —00, 00) we arrive at the equation 

d m—5a(r) 
dr 

\ dr r Po + m c 

where we introduced the notation 



9 _ „gU„M * ' I *( r ) = . < 18 > 



*w= T us) 

\92{r)J 
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for the solution $(r). 

The solutions in the exterior space, r > R, are 




Jm-s+i(kr) 
J m ^s(kr) 

y /po + m c Js- m -i(kr) 
■y/p - rn e J 5 -m(kr) 



for m + 1 — 5 > 

(20) 

for m — 5 < . 



The Jost solution of Eq. ([TJ]) is the solution which behaves for r — > as the free 
solution (pOP . Its asymptotics for r — > oo can be written as 

Hr)- 1 - (f m (k)^ (r) + Uk)^ Hil) (r)) , (21) 

where $#(1,2) are the solutions ( p0|) with the Hankel functions instead of the 
Bessel function. The coefficient f m (k) is the Jost function and f m (k) its complex 
conjugate. 

The ground state energy can be expressed in terms of the Jost function much 
in the same way as in the scalar case |T3]] . However, due to the translational 
invariance in the direction parallel to the flux tube, we have the energy density 



-oo i n ^) 

instead of the general formula (|6]). After carrying out the integration over k 3 we 
arrive at 



fi' s i r(s-i) 2 



2 2v^Fr(s-l/2) 



(«) 



where (n) denotes the remaining quantum numbers in the plane perpendicular 
to the axis of the tube. Now the sum over these (n) can be expressed through 
the Jost function and we get to the relevant order in s, 

£ = C S Y, Uk(k 2 - ml) 1 -* — ]nf m (ik) (23) 

with C s = (1 + s(— 1 + 2 ln(2/i)))/(27r). This representation can be obtained in 
much the same way as done in the scalar case in [pl|. One has to take into account 
the known analytical properties of the Jost function (which differ from that in the 
scalar case). One has to use as an intermediate step a finite quantization volume 
with appropriate boundary conditions (bag conditions work well). Then, in the 
course of tending this volume to infinity, the translational invariant contribution 
from the Minkowski space must be dropped and the remaining finite part, after 
a deformation of the integration contour, just delivers fl2"3]). 



7 



We remark that in the considered problem there are zero modes Jl6| (at k = 0). 
In the sum in Eq. (^) they have to be taken into account. Just in the same way 
as shown in detail in the scalar case in for the bound states they do not show 
up explicitely in representation Eq. (p3|) . 

Here we have taken into accout that both signs of the one particle energies as 
well as both signs of the spin projection give equal contributions to the ground 
state energy thus resulting in a factor of 4 which is included into C s . This 
expression will be used in the calculations below. 



The renormalisation of the ground state energy is defined by Eq. (|12|) . The 
remaining task is to perform the analytical continuation as s — > 0. However, this 
is not immediately possible using representation (|23|) for £ . To continue we use 
the uniform asymptotic expansion In f^(ik) of the logarithm of the Jost function, 
In f m (ik), defined in such a way that the difference 

\nf m (ik)-\nf™(ik) = o(^-) (24) 

is of the order m -4 in the limit m — > oo, k — > oo for *r fixed. Then we split the 
renormalised ground state energy by adding and subtracting \nf^(ik) to get 

gren = £ f + £as ^ 

with the 'finite' part 

1 00 °r 8 
£i = ^ E J dk(k 2 - m 2 e ) — (In f m (ik) - ln/»(i*)) , (26) 



m=— oo 



where it was possible to put s = 0, and the 'asymptotic' part 

d 



£<* = Cs y dA: (k 2 - m 2 c ) l - s — In f™(ik) - £ div . (27) 



m=— oo, 



Here, we included the subtraction of £ dlv according to (|12"D. The continuation 
to s = will yield a finite result because the pole contributions cancel. This 
continuation will be done in section |5|. But because the expression for In / as is 
quite simple, this task can be done analytically. 



The subdivision fl2"5|) of £ is not unique, only the condition (p4]) has to be sat- 
isfied. The inclusion of higher orders into In / as would speed up the convergence 
of the momentum sum in £ f , for instance. But we will use ln/ as in the minimal 
form obeying ([24]). 



S 



4 The uniform asymptotic expansion of the Jost 
function 

The uniform asymptotic expansion of the Jost function can be obtained from the 
Lippmann-Schwinger equation in much the same way as it was done in the scalar 
case |TH]. We rewrite Eq. (JTST) in the form 



d 



d m 



-6a(r) I 1 



■ rn+! * ' $(r) = — ^ " f (r) e A?(r)$(r). 

x d? r Po + m c J r \ i Uy 

The operator on the lhs. can be inverted using the free solutions (|20| ) and we get 
the integral equation 

r 

$(r) = $°(r) + J dr'r'g(r,r') AP(r')$(r') (28) 
o 

with 

g{ry) = (*jr(r)*£«(0 - ®j(r')<S> T H(1) (r)) , (29) 

where $ T means the transposed of Inserting $j ( p0|) for $°(r) in the rhs. of 
Eq. (|28|) , this equation determines just the Jost solution. Using the asymptotic 
expansion of the Hankel functions and comparing with Eq. ([H]) we obtain for 
the Jost function the representation 

oo 

f m (k) = 1 + ^ J dr r $£ C1) (r)A?(r)$(r) . (30) 
o 

Equation (^) can be iterated. It turns out that we need all contributions up to 
the fourth power in AV in order to satisfy condition (p4|). Note that in the scalar 
case the second power had been sufficient. Iterating Eq. ( |2"8"D we obtain 

$(r) = $°(r) 

r 

+ J dr'r'g(r,r')AV(r')$°j(r') (31) 
o 

r r' 

+ J dr' r' J dr"r"g{r, r')AV(r')g(r',r")AV(r")<S>°j(r") 



o 

r r' r" 

+ J dr' r' J dr"r" J dr"'r"'g{r, r')AV{r') 



xg(r' ,r")AV(r")g(r" ,r"')Ar(r'")^ J (r"') 

+0 ({AV)^ 



This expression has to be inserted into Eq. fl3"0|). In fact, we need the logarithm 
of the Jost function. Therefor the appearing expression must be expanded again. 
We write the result as \nf m (k) = J2 n >i m fm HQ where n denotes the power 
of the operator AV. Up to the fourth order we obtain (for several details see 
Appendix A), 

oo 

ln/W(A;) = (^yjdrr^ w (r)AV(r)Qj(r), (32) 
o 

2 oo r 

lnf£\k) = -g) Jdrr J dr' r' $^(r)A?(r)f ffW (r) 

o o 

x^(r')AV(r')$j(r'), (33) 

2 oo r r' 

]nf®(k) = Jdrr J dr' r' J dr" r" (r) AV(r)$ Hm (r) 



x^ T HW (r , )AV(r , )^j{r')^(r'')AV(r")^ J (r''), (34) 

^ oo r r' r" 

lnfg\k) = - Jdrr J dr' r' J dr" r" J dr'" r'" 



(4$£ (1) (r) AV(r)^ HW (r) (r') AV{r')<$>j{r') 





$^ 1) (r")A?(r")f J (r")^(r"')AP(r"')$ J (r'") 
+2^ (1) (r)AP(r)$ H(1 , (r) (r')A?(r')$ H(1) (r') 

^(OATVOMO^KO A?(/")$j(0) > ( 35 ) 

where rearrangings of the integration domains had been made. 

Now, because we are interested in the Jost function for imaginary momentum, 
we turn from the Bessel functions to the corresponding modified ones. Then we 
have to perform the uniform asymptotic expansion of these expressions. 

Before doing this we note that it turned out to be more convenient not to use 
the orbital momentum m as the expansion parameter, but instead 

f m+i for m = 0,1,2,... 

v =\ i f , o ( 36 ) 

[ — m — g tor m — — 1, — 2, . . . 

with z/ = |, |, . . . in both cases. 

Then we need uniform asymptotic expansions of the modified Bessel functions 
for v — > oo, z fixed, of the following type, 



W/2(w)J V 2zy U + t/ i 24 z/ 

-4et 3 -4t 4 + 5et 5 + 5t 6 
+ 



16 z/2 
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-2160 1 4 - 2304 et 5 + 7440 1 6 + 7695 et 7 - 5400 1 8 - 5525 et 9 
+ 5760 v 3 



with e = =Fl for ^ respectively 7^, r\ = y/1 + z 2 + ln(z/(l + yl + z 2 )) and 

t = (1 + z 2 ) 2 which may be derived from the commonly known one (I7j by 
a corresponding reexpansion. Similar expansion for K v -i/2 and I v ~\/2 are also 
used. 

Now we insert this expansion into the logarithm of the Jost function. Then the 
integrations over r'", r" and r' can be carried out successively by the saddle point 
method (only equal arguments in the functions f] (z) yield contributions which do 
not exponentially decrease for v — > oo) as done in |13|]. In doing so it becomes 



apparant that terms up to the fourth power in AV contribute to the asymptotic 
expansion in v up to the order v~ 3 . The relevant saddle point expansion is 
presented in the Appendix B, see (|62|). Finally we collect all contributions up to 
this order and define 



3 3n 00 



/Hr ft 
~ X ^~ n (37) 

n=i,=n n ' V 



with the notation t — (1 + (rk/u) ) 2 . The coefficients turn out to be 

v (aS) 2 v _ (a8) 2 

X 2i2 = \5 2 {a 2 - raa') , X 2>4 = \5 2 (-3a 2 + raa') , X 2j6 = \{a5) 2 

X 3i3 = \5 2 (a 2 -raa' + \r 2 aa" -\5 2 a A ) , (38) 

X 3 , 5 = |5 2 (-fa 2 + 7raa' -r 2 aa" + 65 2 a 4 ) , 

X 3>7 = \5 2 (35a 2 - 5raa' - 5<5 2 a 4 ) , X 3 , 9 = =§5 2 a 2 . 

Here, a means the profile function a(r) in (§). Below, when inserting this ex- 
pansion into £ as , the sum over the orbital momentum must be performed. There 
some contributions cancel, for instance those which are proportional to 6 and S 3 . 
They are not shown in formula (B3). 



5 The asymptotic part of the ground state en- 
ergy 

The asymptotic part of the ground state energy is given by Eq. (p7|) and the 
expression fl3"7] ) for the asymptotic expansion of the logarithm of the Jost function. 
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We rewrite it in the form 



2C„ 



_ 1 3 



dk k 



m. 



n=lj=n 



(39) 



Here, the sum over the orbital momentum m in (B7|) is rewritten as a sum over 
v, Eq. fl36|) . By means of (|63|) it will be replaced by two integrals. The contribu- 
tion resulting from the first integral can be calculated explicitely using formula 
( 64 ) . Together with the explicit expressions for X n ^ (^8|) , after a straightforward 
calculation, it can be seen to cancel exactly £ dlv (for arbitrary profile function 
a(r)). 

So we are left with the contribution resulting from the second integral in 
There the integration over k can be carried out by means of formula (|65j) . Then 
defining 



S nj -(x) 



—i 



dv 



xi J 1 + exp(27rz/) 



(40) 



! + (f) : 



2\ "fa" 



1+ ^ 



2\ A ^2' 



we arrive at 



i 3 3n -:j 

£ as = EE/ s - 

^ n=l i=n ^ 7 



rm e ). 



(41) 



In the functions S nj (a;) 



the analytic continuation to s = has to be per- 
formed. For this reason one has to integrate by parts several times to get rid of 
the singular denominator. The resulting expressions are shown in the Appendix 



C, formulas flBlf ). Now in the integration over r, it is useful to integrate by parts 
which reduces to the substitutions raa' — > —\a l rd r and r 2 aa" — > —a' 2 + \o?rd 2 r. 
We note that after doing this the contribution resulting from n = 2 vanishes, 
already before the integration over r will be carried out. Then the integrations 
over r and v can be interchanged and after rescaling v — > vrm e we get 



S 3 



— 16 r dr 



71 



J [a(r) 2 gi{rm c ) -r 2 a(r)' 2 g 2 {rm c ) + a(r) 4 g 3 {rm c )} , (42) 



with 



9i (x) = / dv Vu 2 ^ 2 fi(u) (i = 1, 2, 3) . 
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The functions fa are displayed in the appendix, Eq. (p7|) . This is the final formula 
for £ as for an arbitrary profile function a(r). 

For the homogeneous magnetic field inside the flux tube, i.e., for the profile 
function (|j), the integration over r can be performed explicitely. After elementary 
calculations we get 



£ £ 



7TR 2 



dz/ 



+ 



dz/ 



Rm a 



3(Rm e ) 2 



35 



(43) 



z/ 3 — Jv 2 



(Rm e ) 2 



3(Rm c 



+ 



(Rm e 



2v 



m e R 



-Af 2 (v)5 2 



v 2 — (Rm e ) 2 



3(Rm c 



x 



8z/ 7 - ^z/ 2 - (Rm c ) 2 (8z/ 6 + Au 4 (Rm e ) 2 + 3u 2 (Rm e ) 4 - 15(Rm e 

105(Rm e ) 6 

(Rm e ) 2 , (v + Jv 2 - (Rm c 
In - — 



+ 



z/- 



(Rm e 



2 2u m e R 

This expression consists of two parts which we write in the form 



£ 3 



2 ei(Rm e ) 4 e 2 (i?m e ) 
o — h o 



(44) 



R 2 R 2 

Here e\ respectively e 2 describe the contributions proportional to the second 
respectively fourth power of the coupling 5 to the background. They are shown 
in Fig. 1. Their behaviour for x — ► can be calculated from (|43| ) and we have 

In x 



ei(x) 
e 2 (x) 



+ 0.1348 + 0(x) 



3tt 

-0.0354 + 0(x) . 



(45) 



The logarithmic contribution is just that which was to be expected from the heat 
kernel expansion (|i~5p. 



6 The 'finite' part of the ground state energy 
and numerical results 

The finite part of the ground state energy is defined by Eq. (|26|) together with the 
asymptotic expansion of the Jost function, Eq. (|3"T|). In general, these quantities 
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-0.04 



-0 . 05 L 

Figure 1: The functions e\ and e2 appearing in the asymptotic part of the ground 
state energy 

can be calculated only numerically. We consider here the case of a homogeneous 
magnetic field inside the tube as given by Eq. ([|). In that case the solutions 
of the field equations are known, they are hypergeometric functions inside and 
Bessel functions outside. As these formulas are in general quite well known resp. 
easy to derive, we give here only the result. The notations are close to that in 
the paper [18]. For positive orbital momentum (m = 0, 1, ... in formula fl5Ep) we 
have 

, _ / kR\ u+1/2 exp(— 5/2) [ kR , f m / {kR) 2 3 \ 



+ v + 2 ) K ^h-^ kR ) iFi rjr> v + 2 ; 6 ) s ' (46) 

and for negative m (m = —1, —2, . . .) 

+ { v + \) ^ + l 2+S (kR) iF, + \ + M! i/ + I; ^ I . (47) 
Thereby the pure Aharonov-Bohm phase is dropped as it does not contribute to 
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Figure 2: The function R 2 5 2 £ i (R) for several values of 5 



The asymptotic part of the Jost function can be obtained explicitely by car- 
rying out the elementary integrations over r in ([J7|). Now, having given all 
ingredients in the integrand of £ f , the remaining task is to perform numerical 
computations for several values of the parameters. For this task it turned out to 
be useful to integrate by parts and to substitute k = y/x/R. Then we have 

oo 

£ f = ^4 E / dx (ln/+(ifc) + ln/-(ifc)-21n/ M (ifc)) .(48) 

This expression can be calculated numerically. The integration over x is quite 
quickly convergent, the sum over v not. So, in order to achive a satisfactory 
precision for the plots, in depedence of the parameters v must be summed up to 
15, for large x even more. 

The general behaviour of £ l as a function of the radius R of the flux tube is 
quite smooth. For R — > it is proportional to R~ 2 . This can be seen analyti- 
cally from Eq. fl48|). For R — > oo it is proportional to R~ 3 which we observed 
numerically. Having in mind that the behaviour for R — > oo is determined by 
the next heat kernel coefficient after a2 we conclude from this that 05/2 is non- 
vanishing. This seems in contradiction with the general results saying that for 
manifolds without boundary half-integer coefficients vanish. But one has to re- 
member, that higher coefficients contain (at least) squares of derivatives of the 
background field which for the presented example leads to undefined expressions. 
Thus, for the higher coefficients the general formulas do not apply and there is 
no contradiction at all. 

In Fig. [|the function £ { {R) is shown multiplied by R 2 5~ 2 as a function of R 
for several values of 5. 

In Fig. |3] the complete ground state energy, £ ren (R), multiplied by R 2 5~ 2 
is shown for several values of 5. In general, this function takes only negative 
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Figure 3: The complete ground state energy, multiplied by R 2 S 2 for several 
values of 5 

values, relatively weakly depending on the flux 5. For small R, the logarithmic 
contribution is dominating. 

The complete energy is the sum of £ class ( ]TB| ) and S ren Q2"5|). In Fig. 3, the 
classical energy would be a straight horizontal line at 2n/a. From this it is 
clear that the complete energy, remaining a monotone decreasing function of the 
radius, deviates only slighly from the classical energy for all values of the radius 
R except for very small ones as mentioned in the end of section 2. 

For large 5, in £ l and £ as the contributions proportional to S 4 dominate, giv- 
ing (at last for < -R-l) £ f resp. £ as large positive resp. negative values. But 
these contributions cancel each other. This was seen in the numerical calcula- 
tions. Also, this corresponds to the procedure of adding and subtracting In / as in 
section 3 which contains terms proportional to 5 4 . However, we did not perform 
a complete investigation of the behaviour for large S. 

7 Conclusions 

In this article we have provided a full analysis of the ground state energy of the 
spinor field in the background of a straight magnetic flux tube of finite radius R. 
The formalism developed applies in principle to any magnetic field with cylindri- 
cal symmetry. Assuming that the Jost function is known or can be determined 
numerically, Eqs. (|42]) and (^) give the final formulas for this case. We have 
applied these formulas to the case of the magnetic field (|2]). The final result con- 
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sists of a very explicit "asymptotic" part, Eq. ([43]), and a part to be determined 
numerically, Eq. fl48|). A detailed numerical analysis shows, that the ground state 
energy turns out be negative, remaining for almost all values of the radius R by 
a factor proportional to the fine structure constant a smaller than the classical 
energy As a result, in the range of applicability of our results, the total energy 
remains positive and, furthermore, does not show a minimum for finite values 
of R given a fixed flux. The magnetic string thus remains unstable also when 
including quantum corrections into the total energy. 

An interesting question is on the general dependence of the ground state en- 
ergy on the specific background chosen here. First we remark, that we are in 
agreement with Q with respect to the ground state energy beeing a small ad- 
dendum to the classical energy (as fas as R is not too small). Secondly, from the 
explicit formulas like (^) for £ as we expect that the dependence on the shape of 
the magnetic field inside the flux tube will be weak. This will be so at last for 
sufficiently smooth background fields. When, however, the background becomes 
singular, this may change. As an example one can consider the ground state 
energy of a massive spinor field with bag boundary conditions on a sphere calcu- 



lated in |2(J showing as function of the radius even changes of the sign. Another 
consequence of our calculations is that it seems impossible to shrink the radius 
of the string to zero bcause of the logarithmic singularity (|i~5|) appearing in that 
case. In view of this, it would be interesting to reconsider earlier investigations 
in the background of the infinitely thin string whereby we admit that it might 
well happen that this singularity can be absorbed into some counterterm. 

Even if we have found a negative answer within the class of examples (§), the 
results presented can be a starting point to consider further the question if inho- 
mogeneous magnetic fields can minimize the energy for fixed flux. Furthermore, it 
seems possible to include other aspects as for example external electric fields and 
the anomalous magnetic moment. Also, the techniques developed here, are suited 
for the calculation of the fermionic contribution to the vacuum polarisation in 
the background of the Nielsen-Olesen vortex, the Z-string or in a chromomagnetic 
background. 

It might be interesting to compare the ground state energy found here, which 
is complely negative, with the effective potential in a homogeneous background 
field which was mentioned in the introduction. For small fields it takes positive 
values which are important for the electroweak symmetry restoration. So our 
result suggests that this might be altered by a finite spatial extension of the field. 
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Appendix A: Perturbation theory for the loga- 
rithm of the Jost function 

In this appendix we will derive the expansion for the logarithm of the Jost func- 
tion, Eq. fl30|). The first step is to use Eq. (|3l| ) in Eq. (|30"1). The Jost function 
itself up to the fourth power in the perturbation (A(P) reads 

f m {k) =:l + x 1 +x 2 + x 3 + x 4 + ((AP) 5 ) , (49) 
with the definitions 

7T 



2z, 



) J™ d n n <$> T m (ri)AP(n)$j(n), (50) 

(yj- \ roc rri 
- J dr ^ r i J Q ^2r 2 ^i(ri)A?(r 1 ) 9 (r 1 ,r 2 )AP(r 2 )$ J (r 2 ),(51) 

x 3 = (— J J dr 1 r 1 dr 2 r 2 dr 3 r 3 <b H i{r l )iSP{ri) 

g{ ri ,r 2 )AV{r 2 )g{r 2 , r 3 )AV(r 3 )$j(r 3 ), (52) 
Xi = (t^tJ J dr x r\ dr 2 r 2 dr 3 r 3 dr A r 4 (r i )A'P(r 1 ) 

^(n, r 2 )AP(r 2 )#(r 2 , r 3 ) AV(r 3 )g(r 3 , r 4 )AV(r 4 )$j(r 4 ).(53) 
We will need the combinations 

ln/ m (*) = In /£>(*) + lnf^(k) + In /,<?(*) + In /£>(*) + O ((A(P) 5 ) , (54) 
with 

ln/W = Sl> (55) 
ln/i 2 ) = xa-ixf, (56) 

hi/i 3) = -x\ - Xl x 2 + x 3 , (57) 

In = --x\ + x?x 2 - -x\ - x x x 3 + x 4 . (58) 

Let us consider ln/ m (fc) order by order. The first order Inf^ik) is already 
given by its definition, Eqs. ([55]) and fl5"0"|). For the calculation of \nfffl(k) some 
manipulations are needed. The main trick, also for the calculation of the higher 
orders, is the rearrangement of integration domains. At the beginning we will 
give details, later on only an idea of the single steps is given. 
Using Eq. (|29|) one obtains 

\nf£\k) = (^) 2 {/°°^m ^ ri ^ 2 r 2 $^ 1 (r 1 )AP(r 1 )$ J (r 1 )<l>l, 1 (r 2 )AP(r 2 )$ J (r 2 ) 



/ drm / rfr 2 r 2 $^ 1 (r 1 )AP(r 1 )$Hi(ri)$J(r 2 )AP(r 2 )<l»j(r 2 ) 
Jo Jo 

/ driri^rOAPtri^jtn) 
Jo 



The first and third term combine to give 

1 /7T 



2 \2iJ Jo 

n 



poo 

/ rfriri$^ 1 (ri)AP(ri)$ J (ri) x 
J o 

/■ri poo 

/ dr 2 r 2 ^(r 2 )A?(r 2 )$j(r 2 )- (ir 2 $^ 1 (r 2 )AP(r 2 )$ J (r 2 ) 

JO Jri 



Next the integration domains may be rearranged, 

fOO rT2 



poo poo poo pT2 

I dr x r\ \ dr 2 = \ dr 2 I dr\ 

JO Jri Jo Jo 



and changing finally the name of the variable, T\ <-> r 2 , one arrives at Eq. (33). 

When calculating the higher orders it is extremely helpful to systematically 
use the lower orders already obtained. So for the next order we start with 

ln/i 3) W = *s - l^f^) 3 ~ On/ffXln/W), 

where 

^3 = 373,1 + ^3,2 + ^3,3 + ^3,4 

consists of the pieces 

' 7T \ 3 f°° f 7 ' 1 f T2 rp 

—J dnn dr 2 r 2 dr 3 r 3 $ H i (n) A?(ri)$j(ri) 
$^ (r 2 ) AP(r 2 )$ j(r 2 )$£i (r 3 ) AP(r 3 )$ j(r 3 ) 

7j- \ 3 poo rr± rr2 



X3A 



X 3 2 

\2i 



/•CO /*T2 

/ dnn / dr 2 r 2 / dr 3 r 3 (n) AP(ri)$j(n) 

JO JO JO 



$^ (r 2 )A?(r 2 )$ H , (r 2 )$j (r 3 ) AP(r 3 )$ j(r 3 ) 

/ yj- \ 3 <>0O /-n /T2 

x 3j3 = - — / dnn dr 2 r 2 dr 3 r 3 (n) A?(ri)$/fi (ri) 
Viz/ Jo jo jo 

$ J (r 2 ) AP (r 2 ) $ j (r 2 ) ^ (r 3 ) AT (r 3 ) $ j (r 3 ) 

/ 7T \ 3 f°° f'2 rp 

x 3A = I— J drin J q dr 2 r 2 y dr 3 r 3 $ H i (n) AP(ri)$ H i (n) 
$J(r 2 ) AP(r 2 )% (r 2 )$J(r 3 ) AP(r 3 )$j(r 3 ) . 

Cancellations occur due to the identity 



oo rri rr n \ 

dnn / dr 2 r 2 ... / dr n+1 f(r 1 )...f(r n+1 ) 



lo Jo Jo [n + ly. 

which can be proven by induction. It shows 

Z3,i~(rn/«) 3 = 0. 
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oo n n+l 







drf(r) 



(59) 



To manipulate the contribution (In /^)(ln f^}) integrals are spitted according to 

oo rV2 roc 

dr 3 = dr 3 + dr 3 (60) 
o Jo Jr 2 



and identities of the kind 



pr poo pr pT2 poo pr 

/ dr\ I dr 2 = / dr 2 / rfr v + / dr 2 \ dr-y (61) 

JO Jri Jo Jo Jr Jo 

are used. One arrives at 

- fafm)0nfm) = _X 3,2 - X 3>3 + £ 3 , 4 

ending up with Eq. (|35|) . 

Finally, the last order we will need can be written as 

ln/i 4) W = ^-l(ln£\k)) 2 -l(\nf^(kf(ln£\k)) 
"(ln/i 1) (A;))(ln/i 3 )(A;))-l(ln/i 1 )) 4 . 

The contribution 

8 

X4 = X 4) j 

8=1 



consists of 



poo pri pv2 pT3 

x 4A = / dr l r 1 / dr 2 r 2 / dr 3 r 3 / dr 4 r 4 $Ji(ri)AP(ri)$j(ri) 

JO jO JO JO 

$£1 (r 2 )AP(r 2 )$ 7 (r 2 )$^ (r 3 ) AP(r 3 )$j(r 3 )$^i (r 4 ) AP(r 4 )$j(r 4 ) 
/•oo r»i rr-2. rr-s 

x 4i2 = - dr 1 r 1 dr 2 r 2 / dr 3 r 3 / dr 4 r 4 $ H i(ri)AP(ri)$j(ri) 
Jo Jo Jo Jo 

$£1 (r 2 )AP(r 2 )$j(r 2 )^ (r 3 ) AP(r 3 )$ H i (r 3 )$J(r 4 )A?(r 4 )$ J (r 4 ) 

/•OO rri /"T2 /T3 

x 4j3 = - dr x ri dr 2 r 2 / dr 3 r 3 / dr 4 r 4 $^i(ri)A7 ? (r 1 )$ J (ri) 

JO JO JO JO 

$£1 (r 2 )A?(r 2 )$ fl i (r 2 )$ J(r 3 ) AP(r 3 )$ j(r 3 )$^ (r 4 ) AP(r 4 )$ j(r 4 ) 

/•oo rr 1 rr2 rr 3 

x 4j4 = / drirx I dr 2 r 2 I dr 3 r 3 / dr 4 r 4 § H1 (r 1 )AV(r 1 )§ J (r 1 ) 
Jo Jo Jo Jo 

(r 2 )A?(r 2 )$^ (r 2 )$J(r 3 ) AP(r 3 )$ H i (r 3 )$J(r 4 ) AP(r 4 )$ j(r 4 ) 

/•OO /-n /T2 rrz 

x 4j5 = - drxn dr 2 r 2 / dr 3 r 3 / dr 4 r 4 $^i(ri)A'P(r 1 )$ H i(r 1 ) 
Jo jo jo jo 

$J(r 2 )AP(r 2 )$ J (r 2 )$^ 1 (r 3 )AP(r 3 )$j(r 3 )$^ 1 (r 4 )AP(r 4 )<l>j(r 4 ) 
/■oo m rra /T3 

x 46 = / driri / dr 2 r 2 / dr 3 r 3 / dr 4 r 4 $ H i(ri)A?(ri)$ ff i(ri) 
Jo jo jo jo 

$J(r 2 )AP(r 2 )$ J (r 2 )$^ 1 (r 3 )AP(r 3 )$ H i(r 3 )<l>J(r 4 )AP(r 4 )$ J (r 4 ) 
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poo m rr2 r-rz 

x 4J = / dnn / dr 2 r 2 / dr 3 r 3 / dr 4 r 4 ^ H i(r 1 )AV(r 1 )^ H i(r 1 ) 
Jo Jo Jo Jo 

$J(r 2 ) AP(r 2 )$ H i (ra^JfcO AP(r 3 )$j(r 3 )$£i (r 4 )AP(r 4 )$j(r 4 ) 
/•oo m /T2 m 

x 4j8 = - dr 1 r 1 dr 2 r 2 / dr 3 r 3 / rfr 4 r 4 $ Jf i(ri)AP(ri)$ H i(ri) 

JO JO JO JO 

$ J(r 2 ) AP(r 2 )$ H i (r 2 )$J(r 3 ) AP(r 3 )$^i (r 3 )$J(r 4 )A?(r 4 )$j(r 4 ) . 
Eq. (|59|) shows 

^4.1 -^On = 0. 

With the help of rearrangements as (|6"0|), fl6"l|), and similar ones, it can be shown 
that 

- (In /«(A;))(ln /£>(*)) = -2x 4 , 7 + 4x 4 , 8 - 2x 4 , 4 , 

-^(^^WfQn/m^)) = -^4,5 + ^4,7 - ^4,3 ~ ^4,8 + ^4,4 ~ 2x 4 , 2 , 

and 

roc rn r-r 2 rr 3 

-2 / driri / <ir 2 r 2 / dr 3 r 3 / dr 4 r 4 (r x ) AP(ri)$ H i (r x ) 
Jo jo jo jo 

(r 2 )AP(r 2 )$„, (r 2 )$J(r 3 ) AP(r 3 )$j(r 3 )$J(r 4 ) AP(r 4 )$j(r 4 ). 
Putting all pieces together one arrives at (p5|). 

Appendix B: Saddle point expansion of integrals 



For the derivation of Eqs. fl37|) and (|38|) repeated use of saddle point expansions 
was made. The relevant result is stated in this Appendix. 

For v — > oo one obtains the following asymptotic expansion, 



oo 

dr' ^{r'y^ = e^ (r) £ (62) 

fc=i 



where the needed leading terms of the expansion are 
Mr) 



h 
hi 
h 2 



(p'(r) 



(<p'(r)f (<p>(r)r 

4>"(r) 3(p'(r)(p"(r) 3<p{r){ V "(r)) 2 <f>(r)(f/"(r) 
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Appendix C: Representation of sums as integrals 

Here we display some formulas used in this paper, where sums are replaced by 
integrals. 

The sum over half integer numbers can be represented by two integrals, the 
necessary analytical properties of the function f{y) beeing assumed: 



<= u 

The first part of £ as in section ^| can be calculated using 



00 *, r ^ - ^y-V- = - m2 ~ 2sn r(2 - * ww* + *=*) 



(63) 



Jm dk v n 2 (rm)"- 1 T(i/2) 

The integration over k in formula (|39"D can be done using 



dk{ e - nrv = _ m2 -, r(2- g) r (s + |-i) fe) 



• (64) 



The functions S nj (x), (f4*0|) can be written as 



! + (-)" 



2\ 



4 



(65) 



(66) 



1 + exp(27rz/) 
/ v 



for n = 1, j = 1, 3 and n = 3, j = 3, 5, 7, 9 with 

/l,3(*> 
/3,s(f 



\1 + exp(27rz/) 

fl ! 

y z/ 1 + exp(27ri/) 
1/1/ v 



3 \u V 1 + exp(27rz/) 



1/1/1 



7/ 



15 If yz/ \ 1 + exp(27rz/) 



1 / 1 / 1 / 1 



105 \u \v \v \1 + exp(27rz/) 
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For n = 2 the formulas are slightly more explicit. They read 



E 2 , 2 = _^ln(l + e - 2 - 
x 

E = 71 1 

3vr 1 7r 2 e 2 ™ 

^2,6 



4x1 + e 2nx 2 (1 + e 2 ™) 2 

The functions fij build the ingredients for the function gi(x) in Eq. (f|2]). Ex- 
plicitely we find 

1 1 1 39 35 35 

= -fi,i{x) - -h,3( x ) + 4/3,3(3:) - — / 3>5 (x) + -g-/3,7(a:) - ^ M x ) 



~^xd x (- J/ 3)3 (a;) + ^/ 3 ,s(V) - 3,7(2) 

+ (|/3,3(x) - |/3,5(x)) , 



M*) = I (fsA x ) - Mx)) , (67) 
Hx) = -I (/ 3)3 (x)-6/3 l5 (x) + 5/3,7(2;)). 



References 

[1] W. Dittrich, W. Tsai, K.-H. Zimmermann, Evaluation of the Effective Po- 
tential in Quantum Electrodynamics, Phys. Rev. D 19, 2929 (1979). 

[2] P. Elmfors, K. Enqvist and K. Kainulainen, Strongly first order elec- 
troweak phase transition induced by primordial hypermagnetic fields, |hcp- 
| ph/9806403[ 

[3] M.P. Fry, Fermion determinants in static, inhomogeneous magnetic fields, 
Phys. Rev. D 51, 810 (1995). 

[4] D. Cangemi, E.D'Hoker and G. Dunne, Effective energy for (2+1)- 
dimensional QED with semilocalized static magnetic fields: A solvable model, 
Phys. Rev. D 52, R3163 (1996). 

[5] G. Dunne, and T.M. Hall, An exact QED 3+ i effective action, Phys. Lett., B 
419, 322 (1998). 

[6] E.M. Serebryany, Vacuum Polarization by Magnetic Flux: The Aharonov- 
Bohm Effect, Theor. Math. Phys. 64, 846 (1985) (Teor. Mat. Fiz. 64, 299 
(1985)). 



23 



P. Gornicki, Aharonov-Bohm effect and vacuum polarization, Ann. Phys. 
202, 271 (1990). 

Yu.A. Sitenko and A.Yu. Babanskii, The Casimir- Aharonov-Bohm effect?, 
Mod. Phys. Lett. A 13, 379 (1998). 

H. Li, D. Coker and A. Goldhaber, Self- consistent solutions for the Vacuum 
Currents around a magnetic flux string, Prys. Rev. D 47, 694 (1993). 

I. Brevik and T. Toverud, Elelctromagnetic energy density around a super- 
conducting cosmic string, Claa. and Quantum Grav,, 12, 1229 (1995). 

C.G. Beneventano, M.De Francia and E.M. Santangelo, Dirac fields in the 
background of a magnetic flux string and spectral boundary conditions, [hep 
th/ 98090811 . 



M. Bordag, Vacuum energy in smooth background fields,. J. Phys. A: Math. 
Gen. 28, (1995) 755-765. 

M. Bordag and K. Kirsten, Vacuum energy in a spherically symmetric back- 
ground field,. Phys. Rev. D 53, 5753 (1996). 

T.B. Branson and P.B. Gilkey, The asymptotics of the Laplacian on a mani- 
fold with boundary, Commun. PDE 15, 245 (1990). 

S.K. Blau, M. Visser and A. Wipf, Zeta functions and the Casimir energy, 
Nucl. Phys. B 310, 163 (1988). 

Y. Aharonov and A. Casher, Ground state of a spin- 1/2 charged particle in 
a two-dimensional magnetic field, Phys. Rev., A 19, 2416 (1979). 

M. Abramowitz and LA. Stegun, Handbook of Mathematical Functions, Natl. 
Bur. Stand. Appl. Math. Ser.55), Washington, D.C.: US GPO, New York: 
Dover, reprinted 1972. 

M. Bordag and S.G. Voropaev, Bound states of an electron in the field of 
the magnetic string, Phys. Lett. B 333, 238-244 (1994). 

M. Bordag, K. Kirsten and D. Vassilevich, On the ground state energy for a 
penetrable sphere and for a dielectric ball, [hep-th/9811015| . 

E. Elizalde, M. Bordag and K. Kirsten, Casimir energy for a massive 
fermionic quantum field with a spherical boundary, J. Phys. A: Math. Gen. 
31, 1743 (1998), |Eep^th/9707083 . 



24 



